Appearance of topological phases in superconducting nanocircuits 
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We construct non-Abelian geometric transformations in superconducting nanocircuits, which re- 
semble in properties the Aharonov-Bohm phase for an electron transported around a magnetic flux 
line. The effective magnetic fields can be strongly localized, and the path is traversed in the region 
where the energy separation between the states involved is at maximum, so that the adiabaticity 
condition is weakened. In particular, we present a scheme of topological charge pumping. 
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When an electron is transported in magnetic field 
around a closed loop, it acquires a phase equal to the 
magnetic flux through the surface spanned by the elec- 
tron path. This phenomenon has been known as the 
Aharonov-Bohm effect for nearly half of a century [1] . In 
the original paper the magnetic field forms a flux line, 
and electrons move in the region of zero field. Berry 
[2] considered this phenomenon to be an early example 
of the geometric phases, which he described in a more 
general context. For each quantum system undergoing 
adiabatic cyclic evolution of its parameters we can find 
phase shifts acquired by its energy eigenstates. Apart 
from the dynamical factor, there is a contribution which 
depends only on the geometry of the path traversed in 
the parameter space (the parameters are time dependent, 
as they are varied throughout the process, but the ex- 
plicit time dependence does not enter the expression for 
the Berry phase, which makes it geometric in nature). 
Berry gives also a formula for this phase in the form of 
an integral of an effective magnetic field (we will refer to 
this field as the "Berry field" ) over a surface spanned by 
the path. This makes the similarity between geometric 
phases in arbitrary quantum system, and the Aharonov- 
Bohm scenario even closer: one can think of the Berry 
field as analogous to the real magnetic field; the param- 
eter space corresponds then to the real position space. 
The Aharonov-Bohm effect in the original setting is, how- 
ever, easier to observe than the Berry phase in general. 
The first reason is that in the former case there is no 
adiabaticity condition constraining the electron velocity, 
while in the latter the rate of the parameters' variation 
should be much smaller than the inverse energy differ- 
ence between the states involved. The second is that for 
strongly localized field any variation in the path of the 
electron does not affect the phase at all as long as the 
path encloses the flux line (this phase is thus topologi- 
cal). For the Berry phase, depending on the system, we 
encounter effective fields usually smoothly varying with 
the parameters, or even uniform, particularly for a spin- 
1/2 system, the phase is the flux of a monopole, or, in 
other words, the area spanned by the traversed loop on 
the unit sphere (see e.g. [2]). Fluctuations of externally 
controlled parameters first of all lead to dephasing of dy- 
namical origin, but also smear the path, which affects the 



visibility of the effect even further. 

Here we show that geometric phases can appear in 
quantum systems which are robust in the sense ex- 
plained above: the effective magnetic field obtained from 
the Berry formula (the Berry field) is not necessarily a 
monopole (uniform in space) field, but can be strongly 
localized near the points of the lowest energy spacing be- 
tween the states involved, and thus the geometric phases 
generated in the process can be topological. The field 
is suppressed in the region where the energy spacing is 
at maximum. This region is thus the most robust, as 
the fluctuations of the parameters do not change the ge- 
ometric phases. Since the energy spacing is large in this 
region, the adiabaticity condition is weakened, and the 
path traversed relatively fast should give strong effects. 

We find here the simplest, two-dimensional non- 
Abelian phases in the system considered in Ref. [3], but 
for a different range of the system parameters. This 
changes the picture substantially, as the topological prop- 
erties of the phases become now evident (geometric 
phases in superconducting nanocircuits have been stud- 
ied also in [4, 5]; in the earlier settings, however, the 
phases have not been proven to be topological). In our 
system the operational subspace is a twofold degenerate 
ground state, so that the dynamical contribution is just 
an overall phase factor, and the actual transformations 
are of purely geometric origin. Since the fields are lo- 
calized, these transformations are topological in nature. 
Furthermore, we find a way to control the localization of 
the Berry field, as well as of the phase acquired while go- 
ing around the region of enhanced field. The position of 
the peaks can be found from the parameter dependence 
of the spectrum; the degree of the localization reflects 
the deviation of the imaginary part of the Hamiltonian 
from zero. In the limit of completely real Hamiltonian, 
the Berry field is zero except for singular points, and the 
geometric phase can have only discreet values {— tt, 0, 7r}. 

The system we consider (see Fig. 1) consists of 
two charge qubits [6] coupled to each other via a 
dc SQUID (superconducting quantum interference 
device). The electrostatic energy of the system, 
including all the capacitances in the system, de- 
pends on the number of the Cooper-pairs on each 
island [denoted by (ni,n%)\, and equals E n \,n2 = 
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FIG. 1: (a) The system composed of two "charge qubits," 
where the state is encoded in the number of Cooper pairs on 
the superconducting islands (black nodes), i.e. m, n.2- The 
voltages and magnetic fluxes can be controlled in a cyclic way 
in order to perform non-Abelian geometric transformations, 
(b) A simplified setup for topological charge pumping. 
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with E cl = 4e 2 C S2 /2(C S iC S2 - C? n ), and sim- 
ilarly for E c2 . Here Csi(2) is the sum of all 
capacitances connected to the first (second) is- 
land. The intcr-qubit electrostatic coupling E m = 
4e 2 C m /(CsiC S 2 - Cf n ). The Josephson coupling J (1) = 
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long to two decoupled subspaces of the Hilbcrt space. 
By varying adiabatically two parameters, e.g. n g2 and 
$2 around a loop, we generate a phase shift between 
those states [exp(z(/>i<7 z /2) in the basis {|0), |1)}] equal to 
the difference of the Berry phases acquired by the states. 
The remaining independent parameter n g \ is adjusted 
accordingly during the variation, so that the degeneracy 
in the ground state is maintained. Similarly for parame- 
ters specified by J\ = J 2 = J and E 01 = E w , the states 
1+) = (10) + |l»A/2 and I-) = (|0) - \l))/y/2 belong to 
decoupled subspaces, and the geometric phase shift be- 
tween them again equals to the difference of their Berry 
phases. In the latter case, however, the corresponding 
transformation is exp(i4> 2 <J x /2), so that the two trans- 
formations generate a non-Abelian group of (all) unitary 
transformations within the subspace {|0), |1)}. 

Let us now recall the original Berry formula for the 
geometric phase. The Hamiltonian is smoothly varied 
in time with its parameters, H[R(t)], and for each point 
R(t) of the parameter space its eigenstates are \n[R(t)]) 
with the energies E n [R(t)). The geometric phase corre- 
sponding to the path C in the parameter space, acquired 
by the state \n), equals then 



|. Here the 

couplings Jd, and J u correspond to the lower and upper 
junctions in the dc SQUIDs, which connect the qubits 
to their reservoirs. The upper junctions J« , and J„ 2 ' 
are here replaced with another dc SQUID loops, which 
enables control over the symmetry of the big dc SQUIDs 
(the meaning of this technique is explained in more 
details in Rcf. [3]). The important fact is that the 
Josephson coupling is real for J u = J^, and complex 

otherwise. Finally J rn = Jm cos 7r<l> m is the Josephson 
energy of the middle SQUID. The magnetic fluxes are 
here normalized to 4>o = hc/2e, the superconducting 
flux quantum. The Hamiltonian of the system reads 

H = E niin2 \nx,n 2 }{nx,n 2 \ - \ [Ji\n\,n 2 ) {n\ + l,n 2 | 

+ J 2 K,rt 2 )(ni,rt2 + M + Jm\n>i + l,n 2 )(rii,n 2 + 1| 
+h.c). (I) 

For the Josephson couplings switched off {J\ = J 2 = 
Jm = 0), the energy eigenstates are the charge states 
with their corresponding electrostatic energies. For the 
particular point n g \ = n g2 = 1/2, the states |0) = 1 01) 
and |1) = 1 10) span twofold degenerate ground state. We 
construct the geometric transformations by choosing an 
appropriate closed path in the parameter space, which 
begins (and ends) at this point (and we refer to it as to 
the starting point). If the gate voltages n g \ and n g2 are 
close to 1/2, only four states are relevant (ni )2 = 0,1), 
and the calculation of the geometric transformations cor- 
responding to given paths is straightforward. 
In particular [9], if we keep the couplings J\ = J m = 0, 
while J 2 complex and finite, the states |0) and |1) be- 
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where the integral is evaluated over the surface spanned 
by C of the effective magnetic (Berry) field 
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The energy difference in the denominator makes it pos- 
sible to predict the regions of suppressed and enhanced 
field. Let us consider the spectral characteristics of our 
system. In Fig. 2(a) we see exemplary diagram of the 
ground-state charging energies. If we move along the bor- 
ders between the cells, the ground state remains twofold 
degenerate. In Fig. 2(b) the spectrum corresponding to 
the thick line in (a) is plotted as the function of n g2 . If we 
now want to perform a phase shift between the states |0) 
and |1), we switch the coupling J 2 to a finite value. The 
coupling will strongly mix the states near the points of 
crossing with the first excited state, as shown in the right 
inset of Fig. 2(b). However, even for finite coupling, the 
minimum in the energy separation between the ground 
state, and the first excited state will correspond to n g2 of 
the triple points in Fig. 2 (a). For these points, the spac- 
ing further depends on \J 2 \, which is tuned by the flux 
<& 2 and has a minimum value for <I> 2 = 1/2 + k and maxi- 
mum for $ 2 = k, k being an integer. The same discussion 
applies to the second transformation, exp(i(j) 2 a x /2). As- 
suming for simplicity that the design is symmetric (char- 
acteristics of the charge qubits are identical), we expect 
minimum spacing for $ = <I>i = <I> 2 = 1/2 + k, and 
n g = n g \ = n g2 corresponding to the triple points ($ m 
is in this case used to maintain the degeneracy in the 
ground state). On the other hand, the field should be 



3 



2 




-1 -0.5 0.5 1 1.5 2 

n g2 




"°' 5 -l -0.5 0.5 1 1.5 2 

ll g2 



FIG. 2: (a) Ground state charging diagram and (b) elec- 
trostatic energies corresponding to the thick line in (a). 
The black node denotes hereafter the starting point. Inset: 
nonzero Josephson coupling mixes the states so that there is 
an avoided level crossing near the "triple points" in (a). 



strongly suppressed (and it is, as shown below) in the re- 
gion of maximum spacing between the ground state and 
the first excited state. Since the aforementioned trans- 
formations are performed by varying both, voltages and 
fluxes, which tune the electrostatic, and the Josephson 
energies, the optimal system should have the character- 
istic Josephson, and electrostatic energies are of compa- 
rable magnitude. [In Ref.[3] the transformations are con- 
structed in the charge regime (electrostatic energies are 
much bigger than the Josephson energies) and the range 
of parameters considered there lies between the minima 
of energy spacing, so that the Berry field seems to be 
smooth and weak.] 

In Fig. 3(a) we see the difference in energy between the 
first excited and the ground state (for J m = Ji = 0, 
E cl = E c2 = E m , 4 2) /E m = 1/2, and J (2) / E m = 2/5), 
as the function of n g2 and $2- The gate voltage n g \ is 
adjusted so that the ground state remains degenerate. 
In Fig. 3(b) we see the Berry field for the same range of 
parameters. The peak position corresponds indeed to 
the minimum in the energy separation. The most robust 
paths should be placed far from the peak, and at maxi- 
mum energy separation (so that the adiabaticity condi- 
tion is weak). In our case the optimal path is the border 
of the rectangle {n g2 ,§ 2 ) € [0,0.5] x [0,1] = S [actually 
the peaks (and dips) form a lattice in the n g2 -& 2 plane, 
so that the most robust paths are borders of all Wigner- 
Seitz primitive cells of the lattice]. One turn along the 




FIG. 3: (a) Parameter-dependent spectrum, and (b) effective 
magnetic field. The minimum in the energy spacing corre- 
sponds to the maximum in the field strength. Variation of the 
parameters around the peak results in a phase shift between 
|0) and |1). (c), (d) Similar characteristics corresponding to 
the phase shift between |+) and |— ). The contours denote the 
most robust paths. 

path produces a phase shift cxp(it/>i<7 z /2), where 

4>i= \ dn g2 d$ 2 V(n g2 ,<!> 2 ). (4) 
Js 

Similarly, for J\ = J 2 = J, and n g = n g \ = n g2 (for 
simplicity we assume that the qubits are identical) we 
vary n g and <E> = <&i = $2 to perform the transformation 
cxp(i02Ca;/2) ($ m is tuned to maintain the degeneracy 
in the ground state), where 

cf> 2 = I dn g d<S>V{n gi $>). (5) 
Js 

The field V(%,$), for J ( d 1)m /E m = 1/2 and 

Ju hi2) /E m = 2/5, is shown in Fig.3(d), and the spec- 
tral characteristics for the same range of parameters in 
Fig. 3(c). The optimal path is here again the border of 
the rectangle (n g , $) <E [0, 0.5] x [0, 1] = S or of any other 
Wigner-Seitz cell of the lattice. 

The phase shift acquired during one cycle can be 
controlled by the strength of the Josephson couplings 
Ji 1),(2) , which are tuned by varying the fluxes through 
the small SQUID loops. Let us illustrate this effect for 
the phase shift 4>\. In Fig. 4(a) we see the phase acquired 
during one cycle as the function of J u /Jd- The phase 
shift increases as the asymmetry of the right SQUID is 
reduced, and tends to the value 7r (and for J u / Jd = 1 is 
exactly 7r as explained below). In Fig. 4(b) we see that to- 
gether with reduction of the asymmetry of the SQUID, 
the width of the peaks tends to zero. Thus the effect 
seems to be amplified together with enhanced localization 
of the field. Simple reasoning shows that for symmetric 
SQUIDs (J u = Jd) the effective field is zero, except for 
singular points. Moreover, loops enclosing an odd num- 
ber of these singularities correspond to the phase shifts 
<f>i = ±7r and 4> 2 = ±7r. Indeed, for symmetric SQUIDs 
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FIG. 4: Dependence of the phase acquired during one turn 
around the region of enhanced field (a) and localization of the 
field (b) as functions of J u /Jd- 



the Hamiltonian in Eq.(l) is real. The Berry field is, 
as seen from Eq.(3) zero, unless there is a crossing in 
the energy difference between the ground state and the 
first excited state. Such singularities are points of zero 
Josephson couplings ($i = $ 2 = &m — 1/2) and can be 
easily identified as the triple points in the charging dia- 
gram (see Fig. 2). Depending on the transformation we 
want to perform, one of the degenerate lowest states (|0) 
or |1) and |+) or |— )) is for each singular point decoupled 
from the first excited state, and the Berry phase will be 
zero for this state, even if the path encloses the singu- 
larity. The complementary state, however, is coupled [in 
the sense that the numerator in Eq. (3) is nonzero] to the 
excited state, and the phase shift will be nontrivial. As 
only two states are involved in the transformation, we 
can now use the spin-1/2 picture. The real Hamiltonian 
H = —B ■ cr/2 is parametrized by the B- field which has 
only x and z nonvanishing components, so that each path 
lies in the plane containing degeneracy. Cyclic evolution 
of the parameters corresponds to solid angles either 0, 
when the loop does not enclose the degeneracy, or ±27r, 
when it does. The corresponding Berry phase is thus ei- 
ther or ±7r (see also [2]). Applying the result to our 
system we see that the phases 4>i and </> 2 can have only 
values or ±7r, so that the resulting non-trivial transfor- 
mations within the subspace {|0), |T)} are ±<j z ,±(t x . 

The latter transformation is a topological charge 
pumping. This can be the easiest effect to observe ex- 
perimentally as it requires a simplified setup [see Fig. 
1(b)], with symmetric SQUID loops only, a single volt- 
age source, and a single current line generating equal 
control fluxes <£>i and $2 (assuming identical qubits). 
Let us estimate the rate of the dynamical dephasing for 



this scheme. Fluctuations of the fluxes $1 = $2, and 
of the gate voltage n g do not violate the constraints 
that we put on the parameters. Thus the most im- 
portant source of errors is in this case the fluctuat- 
ing flux <j> m . Since we are interested in the degener- 
ate subspace only, we apply here the formula for the 
dephasing rate in the environment-dominated regime, 
w 2-KaksT /h [8]. Here a is the constant describ- 
ing the strength of the system-environment coupling. At 
the starting point, where the dependence of J m on the 
flux is strongest, the dephasing rate will be at maxi- 
mum, and we take this value as our estimate. At this 
point a = Rk{MJ^tt/^q) 2 /4:Rj, where the resistance 
of the external circuit inducing the flux, Ri, is typi- 
cally of the order of 100O. Rk ~ 25.8kf2 is the quan- 
tum resistance, and M ~ 0.1 — O.OlnH is the mutual 
inductance between the SQUID and the control circuit. 
The dephasing rate should be compared to the mean en- 
ergy separation between the excited and the ground state 
uj e ^ g , which constraints the time in which the loop is tra- 
versed. Then for ss IK and at T = 20mK we obtain 
T^/uje^g w 10 -5 — 10~ 3 . At optimal design there should 
be thus much space left to both satisfy the adiabatic- 
ity condition and perform many charge pumping cycles 
within the time T^ 1 . 

To summarize, the results presented here prove the 
possibility to perform robust geometric transformations 
in quantum systems by identifying points of localized 
effective magnetic field. The resulting procedure resem- 
bles the original setting of the Aharonov-Bohm effect. 
In this particular system the geometric phases will be 
accompanied by dephasing of dynamical origin, as the 
degeneracy of the ground state is controlled by external 
parameters. However, the weak adiabaticity condition 
on one hand and dominant role of the Berry field from 
the interior of the enclosed region on the other make 
it possible to traverse the paths quickly (and generate 
large phase shifts), before the dephasing destroys the 
visibility of the effect. 
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